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Abstract 

We propose minimal models of one-, two- and three-dimensional micro-swimmers at low Reynolds 
number with a periodic non-reciprocal motion. These swimmers are either "pushers" or "pullers" 
of fluid along the swimming axis, or combination of the two, depending on the history of the 
swimming motion. We show this with a linear three-bead swimmer by analytically evaluating the 
migration speed and the strength of the dipolar flow induced by its swimming motion. It is found 
that the distance traveled per cycle and the dipolar flow can be obtained from an integral over 
the area enclosed by the trajectory of the cycle projected onto a cross-plot of the two distances 
between beads. Two- and three-dimensional model swimmers can tumble by breaking symmetry 
of the swimming motion with respect to the swimming axis, as occurs in the tumbling motion of 
Escherichia coli or Chlamydomonas, which desynchronize the motions of their flagella to reorient 
the swimming direction. We also propose a five-bead model of a "corkscrew swimmer" , i.e. with a 
helical flagellum and a rotary motor attached to the cell body. Our five-bead swimmer is attracted 
to a nearby wall, where it swims clockwise as observed in experiments with bacteria with helical 
flagella. 

PACS numbers: 47.63. Gd, 47.63.mf, 83.10.Rs, 87.17.Jj 
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I. INTRODUCTION 



Micro-swimmers, including motile bacteria and artificial swimmers, swim at low Reynolds 
number by using a periodic non-reciprocal motion, i.e. motion lacking time-inversion sym- 
metry. Since their motion is governed by the Stokes equation, which is linear, a reciprocal 
motion does not induce any time-average net migration of the center-of-mass Jl|. Micro- 
swimmers break reciprocality of the motion, for example, by rotating one or more helical 
flagella attached to the cell body that counter- rotates as a result (e.g. E. coli), or by 
performing "power" and "recovery" strokes or the breaststroke motion with nagella (e.g. 
Chlamydomonas) , or by propagating a wave on a flexible flagellum from the base of the cell 
body to the free end (e.g. sea urchin spermatozoon) These micro- swimmers can be cat- 
egorized into "pullers" , "pushers" or combinations of the two, depending on the swimming 
motion. We here define a puller/pusher as a micro-swimmer that on average pulls/pushes 
fluid to/away from the cell along the swimming axis over one period of its swimming mo- 
tion. Examples of pullers are Chlamydomonas and Chrysomonad, and pushers include the 
sea urchin spermatozoon and E. coli. Some micro- swimmers, such as E. coli and Chlamy- 
domonas, intersperse straight "running motion" with "tumbles" allowing them to reorient 
their movement stochastically towards a more favorable environment or possibly to reduce 
their chances of encountering predators 

The two-bead swimmer has been proposed as a simple model of a micro-swimmer to study 
the collective dynamics of the swimmers and their behavior near a wall 6|, |7| • However, since 
the two-bead swimmer does not account for a periodic non-reciprocal swimming motion, it 
can not generate the time-dependent flow around the swimmer, which is stronger than 

n 

the time-averaged flow [8] and which can induce synchronization of the swimming motions 

of multiple swimmers. Thus, the three-bead swimmer, which is the simplest model that 

r~i 

accounts for a periodic non-reciprocal motion, has also been proposed |9( and used to study 
the interactions of multiple swimmers 10]. However, as we shall describe here, the swimming 
motion in these early studies on the three-bead swimmer is actually a combination of simpler 
puller and pusher swimming motions. 

Here, we propose minimal one-, two- and three-dimensional puller/pusher micro-swimmer 
models, and two- and three-dimensional run-and-tumble swimming motions. By simulations 
using a bead-spring model, we find that the tumbling motion can be easily induced by 
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breaking symmetry of the swimming motion with respect to the swimming axis. 

Finally, we propose a five-bead swimmer as a minimal model for a swimmer with a helical 
flagellum and rotary motor attached to the cell body (i.e. "corkscrew swimmer" This 
five-bead swimmer accounts for the hydrodynamic effects of the rotation of the flagellum 
and the counter-rotation of the cell body as well as the periodic non-reciprocal swimming 
motion. It is found that this swimmer, when placed near a wall, reproduces the behavior of 



E. coli, which is attracted to the wall 
to the wall through the fluid 
tail nor helical motion. 



11] and swims clockwise, when viewed perpendicular 



12), despite the fact that the swimmer does not have a helical 



II. SIMULATION METHOD 

Our minimal model of a one-dimensional swimmer consists of three colinear beads con- 
nected by two bonds whose lengths change in a non-reciprocal manner. This geometry is 
identical to that of previous studies (9), Q, but we induce different cycles (or periodic his- 
tories) of bond-lengths. Examples of non-reciprocal cycles of the two bond lengths, L\ and 
L2, that we impose are depicted in L\-Li space in Fig. [TJ This swimmer can swim at low 
Reynolds number by repeating a non-reciprocal cycle of bond-lengths. 

We confirmed that the bond-length histories in Fig. [TJ induce migration of the swimmer 
by simulations using the bead-spring model. In the simulations, three spherical beads with 



hydrodynamic radius a are connected by two FENE-Fraenkel springs [13J of time-dependent 
equilibrium lengths L\(t) and Lz(t), which follow a configurational history shown in Fig. [TJ 
The FENE-Fraenkel spring force -Fff with equilibrium length L is written as: 

FMr,)=H ( r ; j/L ~L /2 r f (1) 
for (1 - s) < Tij/L < (1 + s), 

where i and j are indexes of beads, is a vector connecting rj (position vector of bead j) 
to Vi, H is the spring constant, and the deformed spring length is restricted to a range 
set by the parameter s, which is 0.01 in our simulations. The motion of each bead of the 
swimmer is computed according to the following equation with the force distribution on 
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beads given by the spring forces 



N 



Vi(t) = ^n ij (r i A'f j , (2) 

or n(t + At) = n(t) + i^Uijirij) ■ /, J At, (3) 

where i>i(i) is the velocity vector of bead i at time t, is the position vector of bead i, 
/j is the force on bead j, N is the total number of beads in the swimmer, At is the time 
increment in the simulation, and "Hy is a hydrodynamic interaction tensor. For Tiij, we use 



the Oseen tensor 
(RPY) tensor 



15 
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with the point-force approximation, or the Rotne-Prager-Yamakawa 



161 ] with more accurate hydrodynamics up to terms of order r- 3 . 
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By superposing the hydrodynamic interactions between beads, this equation with the 
spring forces on beads captures the dynamics of the swimmer. Although both cycles A and 
B in Fig. [TJ induce migration of the center-of-mass at the same speed in the same direction, 
our simulations show that the resulting cycle-averaged flow fields differ. The far-field cycle- 
averaged flow of cycle A corresponds to that of a puller and cycle B to a pusher, and both 
flows decay as r~ 2 . The cycle-averaged flow of the pusher is surprisingly similar to that of 
a detailed bead-spring model of E. coli in a run js] except that the flow of the three-bead 
swimmer does not have angular velocity components around the swimming axis (see Fig. 

By combining the cycle histories of the pusher and the puller, we can create a "square" 
history in L1-L2 space shown in Fig.[TJcycle C, which corresponds to the swimming motion of 
the Najafi-Golestanian swimmer [9j. The cycle-averaged flow of this history is a combination 
of the flows of a puller and a pusher, which decays as r -3 because the dipole components 
of the flow cancel in a cycle that consists of puller and pusher swimming motions. It is also 
possible to generate a net puller or a pusher by setting different ranges for the changes in 



Li and L 2 



10| (i.e. by making the "square" history into a "rectangular" one). 



III. MIGRATION SPEED OF THREE-BEAD SWIMMER 

Here, we evaluate the migration per cycle of a three-bead swimmer whose cycle is repre- 
sented by an arbitrary closed path in L1-L2 space. Once a cycle motion is specified in L1-L2 
space, the corresponding spring-force cycle in /f-/| space is also specified by Eq. [2j where 
f? is the spring force connecting beads i and i + 1, Lj is the time-derivative of L iy and r] is 
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the solvent viscosity: 



// = ft(Lx,L2,Li,L*,a,Ti). 



(4) 



Therefore, we can calculate the migration per cycle L c requiring a cycle time T based on 
the histories of (Li,L 2 ) and (f( , /|) as, 



The migration per cycle does not depend on the cycle time T because of the linearity of the 
Stokes equation, and only depends on the cycle path in L 1 -L 2 space (i.e. L c is an invariant 
of a cycle). For example, if the time spent for a portion of a cycle becomes shorter (e.g. 
x0.5) following the same path, the migration speed becomes faster (e.g. x2). As a result, 
the migration during the portion of the cycle stays the same. 

Now, we first calculate the migration for a small rectangular cycle, which we refer as C T : 
(Li, L 2 ) = (A, B)^(A - 5 U B)^(A -8 U B- 8 2 )->(A, B - 8 2 )^{A, B). The Oseen tensor is 
used for 7Uj in Eq. [2] to calculate the migration (L r ) and we find: 



where g(a, Li, L 2 ) is a rational function (see Appendix A for the derivation), and 8±, 8 2 <C a. 
The contour plot of the function g is shown in Fig. [2j Note that g(a,L\,L 2 ) is always 
positive when the cycle is counter-clockwise (CCW) in L\-L 2 space, and the migration 
along a clockwise (CW) cycle becomes — L r . Since a cycle that does not enclose area in 
Li-L 2 space induces no migration per cycle (i.e. the Scallop Theorem [l|), a cycle can be 
divided into two or more cycles whose migrations sum up to the migrations of the original 
cycle; see Fig. [31 Therefore, an arbitrary cycle can be divided into small rectangular cycles 
whose migrations sum to that of the original cycle, and the migration per arbitrary cycle 
(L c ) can be written as: 



where A is the area within the cycle. The result of L c calculated using Eq. [7] for a CCW 
circular cycle centered at (Li,L 2 ) = (5a, 5a) with radius 0.5a < R c < 3a is plotted in 
Fig. H] along with our simulation results using the Oseen tensor and the RPY tensor. The 
theoretical result from Eq. [7] agrees well with the simulation results. Note that the simulation 
results obtained using the RPY tensor deviates from our theoretical results based on the 




(5) 



L r = g(a,A,B)8 1 8 [ 



"2: 



(6) 




(7) 
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Oseen tensor when R c is large because the RPY tensor, which accounts for more accurate 
hydrodynamics, diverges from the Oseen tensor when the distance between beads becomes 
comparable to the bead radius a. 



IV. DIPOLAR-FLOW STRENGTH OF THREE-BEAD SWIMMER 

By taking advantage of the partitioning of a cycle described in the previous section, we 
here evaluate the cycle-averaged flow field around a swimmer. Once we find the analytical 
expression for the cycle-averaged flow for an arbitrary cycle, we can distinguish between 
puller and pusher swimming from the cycle path in L1-L2 space, since its cycle-averaged 
far-field flow becomes dipolar as following: 



where r is the position vector relative to the position of the center bead of a three-bead 
swimmer, r is the unit vector of r, h is the unit vector pointing in the swimming direction, 
and p represents the strength of dipole, where the flow is that of a pusher when p > 0, that 
of a puller when p < 0, and neither a pusher nor a puller when p = 0. Note that the far-field 
fluid displacement per cycle d(r) = T(v(r)) is an invariant of a cycle. For a small CCW 
rectangular cycle C r , the dipole strength p r can be evaluated analytically in a similar way 
as in the previous section using the Oseen tensor, assuming 5%, 62 <C a and r ;§> a: 



where h(a, L\,L?) is a rational function (see Appendix A for the derivation). Since a cycle 
enclosing no area in L1-L2 space induces no cycle-averaged flow field, we can use the "cycle- 
division" to obtain the dipole strength induced by an arbitrary cycle: 



From the contour plot of the function h(a, L%, L2) shown in Fig. El we find that 1) a cycle 
symmetric across the line L\ = L2 induces no dipolar flow and the perturbation flow decays 
as r -3 , which is the same as a Najafi-Golestanian swimmer, and 2) if a CCW cycle is entirely 
in the region L\ > L 2 , the cycle-averaged flow is that of a puller, and if entirely in L\ < L 2 , 
that of a pusher. The result of p c calculated using Eq. [10] for a CCW circular cycle of radius 




(8) 



p r = h(a, A, B)5i5[ 



'2- 



(9) 




(10) 
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0.5a and centered at (L±, L 2 ) = (LJ, 5a) where 2a < L\ < 20a is plotted in Fig. [6]along with 
our simulation results using the Oseen tensor and the RPY tensor. The theoretical result 
agrees well with the simulation results except that the simulation results with the RPY 
tensor diverges from the theoretical result because of its higher order accuracy compared 
to the Oseen tensor. Note also that if a CW cycle is entirely in the region L\ > L 2 , the 
cycle-averaged flow is that of a pusher, and if entirely in L\ < L 2 , that of a puller. A 
cycle-direction inversion (CCW -H- CW) reverses the migration direction (+x -H- —x), and 
the sign of the flow dipole (p c <r» — p c ) as shown in Fig. [7J This is consistent with the fact 
that p c (or the function h(a,Li,L 2 )) changes sign upon a reversal of the cycle-direction. 



V. RATE OF ENERGY CONSUMPTION OF THREE-BEAD SWIMMER 



The rate of energy consumption during a cycle can be written as: 

e(s) = fm + f 2 v 2 + f 3 v 3 

= mv 2 - Vl ) + ti{v, - v 2 ) = flU + /|L 2 , (11) 

where s is the distance coordinate along a cycle in L\-L 2 space. Due to the linearity of the 
governing equation (i.e. the Stokes equation), each spring force is proportional to the time 
derivative of the spring length, /? = Cj(s)Lj. Using the chain rule of differentiation, we find: 

e(s) = ci{s)L\ + c 2 (s)L 2 
= ci(s) 



dLi\ 2 ( ds\ 2 

-17 U +C2{s) 



dL 2 \ 2 ( ds\ 2 



ds 



dt 



ds 
dt 



Cl ls 



dL\ \ / \ ( dL/2 
+ c 2 [s) 



ds 



ds 



Therefore, 




ds J \ ds 

By integrating Eq{T3] over a cycle, we obtain the cycle time T, 



T 



ci is) 



dL x 

ds 



+ c 2 (s) 



dL 2 
ds 



e(s) ds. 



(12) 



(13) 
(14) 



(15) 



For a fixed cycle in L1-L2 space, ds/dt and T can be derived from Eq. [T3l and [T5l once e(s) is 
specified. Since in the Stokes flow the rate of energy dissipation is quadratic in the velocity, 
It is easy to show that the most energy efficient cycle motion s e s(t) along a cycle in a cycle 
time T is achieved when the rate of energy consumption is constant through the cycle (see 
Appendix B for proof). 

Assuming a constant rate of energy consumption e, the non-dimensional rate of energy 
consumption e and the non-dimensional energy consumed per cycle E can be defined as: 

eT 2 

e = ^3 , (16) 

~ ET , x 

E = — , 17 

■qa 6 

where E = eT, and these non-dimensional parameters only depend on the cycle path in 
Li-L 2 space. Note that by setting e to be a constant through a cycle, its value can be easily 
calculated from Eq. [15] with given cycle time T. 

VI. MULTI-DIMENSIONAL SWIMMERS 

The swimming direction of the one-dimensional three-bead swimmer is restricted to one 
dimension because of its linear geometry. However, we can easily construct a two-dimensional 
swimmer by using beads that are not colinear (Fig. [1]). Some micro-swimmers tumble and 
reorient the cell body by breaking symmetry of the swimming motion over the swimming 
axis. For example, Chlamydomonas swims straight with synchronous beating of two flagella 



and tumbles with asynchronous beating (4], and E. coli swims straight by rotating multiple 
helical flagella in a coherent rotational direction and tumbles by rotating at least one of the 
flagella in the reverse direction |3] . We confirmed by simulations using the bead-spring model 
with time-dependent equilibrium spring lengths analogous to the one-dimensional swimmer 
that our two-dimensional swimmer also tumbles by breaking symmetry of the swimming 
motion as shown in the cycle history in Fig. [TJ Therefore, by alternating between the runs 
and tumbles, this two-dimensional swimmer can swim in a two-dimensional random-walk 
trajectory. 

We can also create a three-dimensional swimmer by using four or more non-coplanar 
beads, for example, on the vertexes of a tetrahedron or "tetrumbbell" This three- 

dimensional swimmer can swim straight with a periodic non-reciprocal motion analogous to 
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that of the two-dimensional swimmer and tumble by breaking symmetry of the swimming 
motion (see Fig. [I]). It is worth noting that a tetrumbbell with constant equilibrium bond- 
lengths but different spring constants can migrate in a shear flow in the vorticity direction 



at low Reynolds number 



17| . This migration (or "swimming" ) also results from a periodic 



non-reciprocal deformation induced by the shear flow on the tetrumbbell. 

All of these models will be useful for studying the low Reynolds number hydrodynamics of 
micro- swimmers, the swimmer- swimmer interactions, and the collective dynamics of many 
swimmers. To study the multi-swimmer interactions of "corkscrew swimmers" such as E. 
coli, however, we need to include the effect of the angular velocity around the swimmer 
induced by the rotation of the flagellum and the counter-rotation of the cell-body. In Fig. |H1 
we depict a five-bead swimmer that induces a time-average flow field corresponding to either 
a pusher or a puller plus this angular velocity. The five-bead swimmer swims via a periodic 
non-reciprocal motion of beads 2, 3 and 4, just as the three-bead swimmer does, by changing 
the bond-lengths L\ and L 2 . Additionally, torques T\ and T2, and counter-torques —T\ and 
— T 2 are applied as shown in Fig. [8] to mimic the rotation of the flagellum and the counter- 
rotation of the cell-body. Note that, for the total torque to be zero, the magnitudes of 
torques Ti and T 2 must be the same, and also that when beads 2, 3 and 4 are colinear, 
torques — Ti and — T 2 cancel out. Each of these torques is first decomposed into two torques 
perpendicular to the bonds of the swimmer, and then each of these decomposed into forces 
on beads. For example, the force distribution on three connected beads i, j and k induced 
by a torque T tot , which represents either Ti, T 2 , — Ti or — T 2 , is (see Fig. El), 

Fi = ^ totX % (18) 

r ij I J- tot x r ij I 

Fj = -F { - F k , (19) 
F k = ^ Ttot X rij (20) 

r jk \Ttot x r ij\ 

To restrict the deformation of the swimmer, bending potentials are applied: 0(# 2 34, 180°), 
0(6*123, #e) and 4>{0 U5 ,8 e ), where <f)(6 ijk ,6 ) = \k h (cos6 ijk - cos6> ) 2 , 9 ijk is the angle formed 
by beads i, j, and k, 9 is the equilibrium bending angle, k^ is the bending potential constant, 
and we choose 8 e = 160° for our simulations. Although the bending potentials allow small 



deformations of the swimmer 



181 ]. neither the time-averaged flow field around the swimmer 



nor the swimming behavior are affected qualitatively for k^ values as large or larger than 
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that used here. The time for L\ and L2 in a five-bead swimmer to traverse one side of a 
triangle in the configuration space is kept constant at the value to- 

By solving Eq. [2] with the force distribution on beads obtained by adding together all 
forces associated with torques and all spring and bending forces, we obtain the time-averaged 
flow field around a five-bead pusher swimmer, which we find to be very similar to that 
obtained from a detailed bead-spring model (with 60 beads) of E. coli with multiple helical 
flagella in a run 8|. Input parameters with physical units in this simulation are R (maximum 
bond- length) , 77 (solvent viscosity), and T (magnitude of torques). Therefore, we scale length 
and time with R and r = rjR 3 /T, respectively. We choose the simulation parameters a=0.2R, 
5=0.5R, t =2.0r, H=1000T/R, k h =2000T, fl e =160° and At=10- 3 r. 

Simulations were also performed for the five-bead pusher near a wall by employing the 
RPY tensor with the wall effect included jjij and a short-range repulsive potential between 
each bead and the wall with cut-off length 0.2R. We find that the swimmer is attracted to 
the wall and swims clockwise (see Fig. [T0|) as observed in experiments and simulations with 
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2l| , despite the fact that the five-bead 



bacteria with helical flagella(um) (e.g. E. coli) [121 . 
swimmer does not have a helical tail nor helical motion. Therefore, the five-bead swimmer 
captures the experimentally observed behavior of E. coli qualitatively, and will allow us to 
simulate the collective dynamics of micro-swimmers more realistically than before at modest 
cost. 

We also confirmed that a three-dimensional tumbling motion of the five-bead swimmer 
can be induced by temporarily connecting beads 2 and 4 with a FENE-Fraenkel spring of 
constant equilibrium length R while the bond-lengths L\ and L2 change in the same way as 
in a run, in other words, by forming a two-dimensional swimmer in a tumble with beads 2, 
3 and 4. Simulations show that the trajectory of two five-bead pushers in a run, initially 
placed side by side and parallel to each other, agrees qualitatively with that obtained from 
a boundary element simulation of a modeled corkscrew swimmer, showing that two parallel 
swimmers first attract each other, and at the same time their swimming axes rotate in 
opposite directions [22[]. We find that two side- by-side pullers behave similarly. 
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VII. SUMMARY AND FUTURE DIRECTIONS 

We have developed simple micro-swimmer models by extending the Najafi-Golestanian 
three-bead linear swimmer model to allow arbitrary history of bond-lengths (Li,L 2 ), by 
making the beads noncolinear to introduce tumbling, and by adding two additional non- 
collinear end beads to which torque and counter-torque are applied to induce a helical flow 
field around it, thus mimicking the swimming of flagellated bacteria. We have shown that 
both the migration distance and optimum energy consumption per cycle is an invariant of the 
cycle path, and can be computed, respectively, by area and path integrals over the cycle. For 
three-bead swimmers, we have also developed a simple criterion to distinguish puller from 
pusher swimmers, based on the cycle in L\-L 2 space. Finally, we have shown that our five- 
bead corkscrew swimmer spirals towards flat surfaces and interacts hydrodynamically with 
other five-bead swimmers in ways qualitatively similar to that of much more refined models 
of E.coli. These minimal models of micro-swimmers will help in the study of the collective 
dynamics of micro-swimmers with a simplified but qualitatively accurate hydrodynamics. 
These models can also be used to study the synchronization of multiple swimmers if the 
bond- lengths are changed by applying time-constant, equal and opposite forces on the beads 
at the ends of each bond until the bond-length reaches a designated length, as proposed in 
a previous work |23|. 
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APPENDIX A: ANALYTICAL EXPRESSIONS FOR THE MIGRATION AND 
THE DIPOLAR-FLOW STRENGTH INDUCED BY A SMALL RECTANGULAR 
CYCLE 

We analytically evaluate the migration per cycle induced by a small rectangular cycle C r : 
(Li, L 2 ) = (A, B)^(A - 5i, B)-t(A -5 ll B- 5 2 )^-(A, B - 5 2 )-*(A, B). Since the migration is 
independent of how the time is spent for strokes and determined only by the path in L\-L 2 
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space because of the linearity of the governing equation, we assume the first stroke to be 
(Li, L 2 ) = (A — jrSi, B) where < t < Ti. By solving the following simultaneous equations 
derived from Eq. [2] with the Oseen tensor, we obtain (f((t), f 2 (t)), 

o f s f s 9 fs I fs fs _ f s 

£ _ ^ 2 _ v = 1 J2 _|_ z il ' J 2 _|_ J2 _| J2 

67rr]a 4irr]L 1 Anr]L 2 47r?](L 1 + L 2 ) ' 

9 f s f s 9 f s \ fs fs f s 

L 2 = ^3 - W2 = =^ /l + fjj + h + /l + Jl _ (22) 

Qnr]a \kt}L 2 AnrjLi Atitj^Li + L 2 ) 

Then, we can calculate (v i(t), v 2 (t), V3(£)) using (fl(t), f£(t)), 

fs fs fs 

Vl = J i + h lA + ^ (23) 

6%r]a AirrjLi A7ir](L 1 + L 2 ) ' 

f s fs fs fs 

V2 = ll^l. + + (24) 

67r?7a AixriLi Atii]L 2 

fs fs fs _ fs 

v ^ _ J2 _|_ Jl J2 _| Jl /25) 

671770 4nr]L 2 Aiir](Li + L 2 ) 
Finally, we find the center-of-mass migration velocity v m i g (^) = + v 2 + v^)/?> to be 

^mig(t) = ^[-2L 1 L 2 (L 1 + L 2 ){(2L! + L 2 )L\ - (L, + 2L 2 )L 2 2 }+ 
3a{(2L a + L 2 )L? + (2L X + L 2 )L?L 2 + + L 2 )L?L2 - (L x + 2L 2 )L X L\ - {L x + 2L 2 )L A 2 }]/ 
{-AL\L\{L x +L 2 f+AaL x L 2 {L x +L 2 ){^ 

(26) 



24j because 



Note that this solution for the migration velocity is different from that in Ref. 
we do not use approximation a/ Li <C 1, which is used in their calculation. 
We find the migration induced by the first stroke using Eq. [261 

/ v mig {t)dt = — {-2AB(-2A 3 ~2A 2 B+AB 2 +B 3 )+3a(-2A 4 -2A 3 B+A 2 B 2 +AB 3 +B 4 )}/ 
Jo 3 

{-4A 2 B 2 {A+B) 2 +4aAB{A 3 +4A 2 B+4AB 2 +B 3 )+3a 2 {A*-2A 3 B-5A 2 B 2 -2AB 3 +B 4 )}+0{5 2 1 ). 

(27) 

By adding up the migrations induced by the other three strokes in the cycle, which can 
be calculated in the same way, we find the total migration for this rectangular cycle to be 
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L r = g(a, A, £>)£i£ 2 , where Si, 62 <C a and g(a, A, B) is 

g(a, A, B) = ^a{lQA 2 B 2 (A + 5) 2 (A 4 + 2A 3 B + A 2 B 2 + 2AB 3 + 5 4 ) + 

36a 2 (A 2 + 3AB + 5 2 )(A 6 + 3A b B + 3A 4 5 2 + A 3 B 3 + 3A 2 5 4 + 3AB 5 + B 6 )- 

27a 3 (A + 5)(2A 6 + 5A 5 B + 6A 4 5 2 - 5A 3 B 3 + QA 2 B A + + 2B 6 )- 

\2aAB(A + J B)(4A 6 + UA 5 B + 17A 4 5 2 + 10A 3 5 3 + 17A 2 5 4 + 14AB 5 + AB 6 )} / 

{-4A 2 J B 2 (A+ J B) 2 +4aA J B(A+ J B)(A 2 +3A J B+ J B 2 )+3a 2 (A 4 -2A 3 J B-5A 2 J B 2 -2A J B 3 +5 4 )} 2 . 

(28) 

The strength of the dipolar-flow induced by the same cycle C r can be evaluated in a 
similar way by calculating the cycle-averaged flow: 

(v(r)) = \l Y.H^r - nit)) ■ Ut)dt (29) 
1 Jo i=i 

where r is the position vector originated at the position of the center bead of a three-bead 
swimmer, and Vi{t) is obtained by ri{t) = Tj(0) + J * Vi(t)dt. We find the strength of dipolar 
flow to be p r = h(a, A, B)5\52 by comparing the result of Eq. [29]with Eq. [8] under conditions 
Si, S2 <C a and r»a, where h(a, A, B) is 

h(a, A, B) = ^cl\-A + B){-12aAB{A + B){2A + B){A + 2B){A 2 + AB + B 2 ) 2 + 

9a 2 (A 2 + 3AB + 5 2 )(2A 3 + 3A 2 5 + AB 2 + 5 3 )(A 3 + A 2 B + 3AB 2 + 25 3 ) + 
AA 2 B 2 (A + 5) 2 (2A 4 + 7A 3 5 + l\A 2 B 2 + 7AB 3 + 25 4 )}/ 
{-AA 2 B 2 (A+B) 2 +AaAB(A+B)(A 2 +3AB+B 2 )+3a 2 (A*-2A 3 B-5A 2 B 2 -2AB 3 +B 4 )} 2 . 

(30) 

APPENDIX B: MOST ENERGY EFFICIENT CYCLE MOTION OF THREE- 
BEAD SWIMMER 

Here, we show that the most energy efficient cycle motion s e s(t) along a cycle in a cycle 
time T is achieved when the rate of energy consumption is constant through the cycle. We 
show this by evaluating the change of the total energy consumed in a cycle when a portion 
of a cycle-motion is sped up and another is slowed down, while keeping the cycle time T 
constant. Let us change the time spent for two small line segments on a cycle, As p and 
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As q , which is passed through during a small time interval AT, keeping the segment lengths 
constant. 

AT — > (1 + e)AT in the segment As p . (31) 
AT ->■ (1 - e)AT in the segment As Q . (32) 

As a result, we find the following changes: 

As 1 As v- f s 

AT^l+^AT' Vl ~*T+l' ft ^l + l in the S6gment ASp ' (33) 

v .^_^ f!^J*~ in the segment As„. (34) 
AT 1 - e AT 1-e' Jl 1-e 6 9 v ; 

These changes cause the changes in the rate of energy consumption, 

3 

e p = ^2 f* v i ~^ 7] _i_ P n 2 in tne se S ment As p , (35) 
i=i l i + e J 

3 

e,, = ^ ->■ 7- _ g , 2 in the segment As g . (36) 



The energy consumed in these segments becomes, 



e„AT ->■ . 6p xo • (1 + e)AT = , 6p . AT in the segment As„, (37) 
p (1 + e) 2 v 7 (1 + e) 



e 9 AT ->■ 69 • (1 - e)AT = Gq AT in the segment As 9 , (38) 

v l t) \L €) 

where e p and e q are the rate of energy consumption before the change of time interval in 
segments As p and As g , respectively If we assume that the original rate of energy consump- 
tion are constant, or e p = e q = e, the change of the total energy consumed in a cycle after 
the change of time intervals is 

f eAT eAT 1 e 2 
AE = \ — + — L - 2eAT = 2^-eAT > 0. (39) 
^ 1 + e 1 — e J 1 — e 2 

Therefore, if the rate of energy consumption changes in the course of a cycle while keeping T 
constant, the energy consumed in a cycle increases. In other words, the most energy efficient 
cycle motion along a cycle in cycle time T is achieved when the rate of energy consumption 
is constant through the cycle. 
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FIG. 1: Minimal models of one-, two- and three-dimensional micro-swimmers and the cycle his- 
tories. Chrysomonad and sea urchin spermatozoon are examples of one-dimensional puller and 
pusher, respectively. Note that tumbling is induced in the 2D- and 3D-swimmers by holding the 
equilibrium lengths of one (2D) or three (3D) bond-lengths fixed while varying the lengths of the 
other bonds. 
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FIG. 2: (color online). Contour plot of function g(a, L\, L<i) in Eq. [7J 
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FIG. 3: The top three figures illustrate that since in Stokes flow no net migration is produced by 
any bead motions that merely reverse themselves, the migration in any single cycle can be obtained 
as the sum of the migrations in two cycles into which the original cycle is divided. Extending this 
principle, the bottom figure illustrates that the migration in a cycle of arbitrary shape approaches 
that of the sum of migrations in small rectangular cycles into which the original cycle is divided. 
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Rc [a] 

FIG. 4: The migration per cycle for a CCW circular cycle centered at (Li,L2) = (5a, 5a) with 
radius 0.5a < R c < 3a, obtained from Eq. [7] (Theory) , simulations with the Oseen tensor and the 
RPY tensor. 
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FIG. 6: The strength of dipole flow, p c , for a CCW circular cycle with radius 0.5a and centered 
at {L\,Li2) = (LI, 5a) where 2a < L\ < 20a, obtained from Eq. [7] (Theory) , simulations with the 
Oseen tensor and the RPY tensor 
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Swim to -x, puller Swim to -x, pusher 



FIG. 7: The change of the swimming behaviors by an inversion of the cycle-direction or an exchange 
of L\ and Li- 
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Five-bead swimmer 




Torque decomposition 




=Ta+Tb 



FIG. 8: (Top) the configuration of the five-bead model. (Middle) the torque distribution on the 
swimmer. (Bottom) a decomposition of a torque Ttot on three connected beads, which represents 
either Ti, T2, — T\ or — T2 in the middle figure, into two torques T a and T& acting perpendicular 
to adjacent bonds. 
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FIG. 9: (color online). The time-averaged flow field around a three-bead pusher (top), around a 
five-bead pusher (middle) and around a detailed bead-spring model of E. coli in a run 8] (bottom). 
All arrows are unit flow-velocity vectors. Red (or light gray) arrows point into the paper, blue (or 
dark gray) arrows the opposite, and black arrows only have in-plane velocities. Swimming is from 
left to right. 
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FIG. 10: A trajectory of the five-bead swimmer near a wall. The trajectory is projected onto 
the wall to show the radius of the circular motion. The swimmer is initially placed 2R away 
from and parallel to the wall. Simulation parameters are a=0.2R, 5=0. 5R, to=2.0r, H=W00T/ R, 
/c b =2000T, 6» e =160° and At=l(r 3 T. 
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